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Abstract 

It is shown theoretically that the Luttinger liquid phase in quasi-one-dimensional conductors 
, can exist in the presence of impurities in a form of a collection of bounded Luttinger liquids. The 



CO 



conclusion is based upon the observation by Kane and Fisher that a local impurity potential in 
Luttinger liquid acts, at low energies, as an infinite barrier. This leads to a discrete spectrum of 



collective charge and spin density fluctuations, so that interchain hopping can be considered as a 
small parameter at temperatures below the minimum excitation energy of the collective modes. 
The results are compared with recent experimental observation of a Luttinger-liquid-like behavior 
■ in thin NbSe3 and TaS3 wires. 
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Electronic properties of one-dimensional (ID) metals are known to be very different from 
those of ordinary three-dimensional (3D) metals (for a review see Ref.[l, 2, 3]). 3D electron 
gas is well described by Landau's Fermi-liquid picture in which interaction modifies free 
electrons making them quasiparticles that behave in many respects like non-interacting 
electrons. In contrast to the 3D case, in ID electronic systems the Fermi-liquid picture breaks 
down even in the case of the arbitrarily weak interaction. In ID metals, the single-electron 
quasiparticles do not exist, and the only low energy excitations turn out to be charge and spin 
collective modes with the sound-like spectrum. These modes are dynamically independent 
giving rise to a spin-charge separation in ID systems. Furthermore, correlation functions at 
large distances and times decay as a power law with interaction dependent exponents. Such 
a behavior has been given a generic name Luttinger liquid [4] . 

The concept of Luttinger liquid is of great interest in view of its application to real 
physical systems, such as carbon nanotubes and semiconductor heterostructures with a 
confining potential (quantum wires and quantum Hall effect edge states). The case of a 
special interest are quasi- ID conductors, i. e., highly anisotropic 3D conductors with chain- 
like structure. Numerous experimental studies of both organic and inorganic qlD conductors 
at low temperatures typically demonstrate broken-symmetry states, like superconductivity, 
spin- or charge- density wave (CDW) states, and a metallic behavior above the transition 
temperature with non-zero single-particle density of states at Fermi energy. For instance, 
the most studied inorganic quasi-lD metals undergo the Peierls transition from metallic 
state either to a semiconducting CDW state (e. g., blue bronze K0.3M0O3, TaS 3 , (TaSe 4 ) 2 I 
etc.) or to semimetallic CDW state (NbSe 3 ) [5]. Typically, this transitions occur in the 
temperature range 50 - 250 K. From the theoretical point of view, the formation of Luttinger 
liquid in quasi- ID conductors at low enough temperatures is also problematic because of 
the instability towards 3D coupling in the presence of arbitrarily small interchain hopping 
[6, 7, 8, 9, 10]. So the interchain hopping induces a crossover to 3D behavior at low energies, 
while Luttinger liquid behavior can survive only at high enough energy scale where it is not 
affected by 3D coupling. 

In contrast to the interchain hopping, the Coulomb interaction between the electrons at 
different chains does not destroy the Luttinger liquid state, the main difference from the 
ID case being the absence of simple scaling relations between the exponents of the various 
correlation functions [11, 12, 13, 14]. 
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However, in recent experimental studies of temperature and field dependence of conduc- 
tivity of TaS3 and NbSe3 in nanoscale-sized crystals a transition from room-temperature 
metallic behavior to nonmetallic one accompanied by disappearance of the CDW state at 
temperatures below 50 - 100 K was observed [15, 16, 17]. The low temperature non-metallic 
state was characterized by power law dependencies of the conductivity on voltage and tem- 
perature like that expected in Luttinger liquid, or by more strong temperature dependence 
corresponding to the variable-range hopping. Resembling dependencies of conductivity were 
reported although in focused-ion beam processed or doped relatively thick NbSe 3 crystals 
[17]. Hopping conductivity was also detected in dirty quasi- ID conductors KCP and organic 
TCNQ-based metals [18], while pure materials are known to undergo the Peierls transition 
to the CDW state. 

In order to account for such behavior, we study the possibility of impurity-induced sta- 
bilization of a gapless Luttinger liquid state in quasi-lD metals. Impurities in Luttinger 
liquid are known to act as infinite barriers forming the effective boundaries for low energy 
excitations [19, 20, 21]. This leads to a dimensional quantization and, consequently, to a 
minimal excitation energy uo\. As a result, the interchain hopping does not destroy the 
Luttinger liquid phase at temperatures T <C Ui, producing only small perturbations of the 
ID picture. To show this, we consider first the gapless ID Tomonaga-Luttinger (TL) model 
with impurity potential included, and make certain that the system with impurities breaks 
up into a set of independent segments described as bounded Luttinger liquid with discrete 
spectrum. Then we calculate corrections caused by the interchain hopping to thermody- 
namic potential and to the one-particle Green's function, and find that such corrections are 
small at low temperatures. Finally, we discuss modifications introduced by generalization of 
the TL model to the more realistic case of Coulomb potential and compare our results with 
experimental data. 

First of all we start with the TL model ignoring interchain hopping integral t± in order 
to formulate the problem in the zeroth approximation in t±. Electronic operators for right 
(r = +1) and left (r = —1) moving electrons with spin s are given in terms of phase fields 
as (see Ref. [1, 3]) 

Vv, s (z) = lim -==^ s e iA % A r = —[Q - r$ p + s(& a - r$ CT )] (1) 
a-*o ^2na v2 

here phase fields <&„(#) are related to charge (y = p) and spin (y = a) densities, while fields 
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Q v (x) are related to the momentum operators U u = (l/ii)d x Q u canonically conjugate to <fv 
Further, rj r>s are Majorana ("real") Fermionic operators that assure proper anticommutation 
relations between electronic operators with different spin s and chirality r, and the cut off 
length a is assumed to be of the order of interatomic distance. 

We describe the intrachain properties of the system by the standard TL Hamiltonian 
[1, 3] with added impurity backscattering term [2]. In the bozonized form it reads 

+ E VodS (x — xi) cos (V2$ p + 2k F x)cos(V2$>a(x)), (2) 

i 

where v v are velocities of the charge {y = p) and spin {y = a) modes, K v = vpjv v is the 
standard Luttinger liquid parameter describing the strength of the interaction, Vq and d ~ a 
are amplitude and radius of the scattering potential, respectively. 

Kane and Fischer [19] found that the backscattering impurity potential for repulsive 
potential (K p < 1) flows to infinity under scaling. Their arguments were generalized by 
Fabrizio and Gogolin [22] to the case of many impurities. It was shown that the impurity 
potential can be considered as effectively infinite provided that the mean distance, /, between 
impurities satisfies the condition 

' » IT ( 77 ) > (3) 



k F \V J 

where D is the bandwidth. We assume that the impurity potential is of atomic scale, Vq <; D, 
and the interaction between electrons is not too weak, (i. e., K p is not too close to 1). Then 
condition (3) is satisfied for I 3> 1/kp ~ a which is of the order of interatomic distance. So 
the limit of strong impurity potential should be a good approximation in a wide range of 
impurity concentrations. 

Further, 11^, 6^ and <fv must obey the commutation relations (see Ref.fl, 2, 3]) ensuring 
anticommutation of electronic operators (1). Using then the analogy of Eq.(2) with the 
Hamiltonian of an elastic string strongly pinned at impurity sites, we can write down solu- 
tions for the phase operators $ y and 0^ in the region between impurity positions at x = Xi 
and x i+ i as 

n=l » n V2<j j<i k 

00 / 1 

®»( x ) = El 77 — ( bn ~ b n) cosq n x + V , (4) 



where x = x — Xi, U = Xi + i — Xi, q n = nn/k, $j is the modulo 2tt residue of 2/cpXj. Further, 
AN pi = (AN Vl + Nn)/V2, AN ai = (AN Vl - N^)/y/2, and AN Vi (A%) is the number of 
extra electrons with spin up (down) in the region between i-th and (i + l)-th impurities, 
and, finally, Q V i is the phase canonically conjugate to AN V i ([9 U i, AN U i] = i). 

Excitation spectra of the eigenmodes are u v = nu>i tU = v u q n where ui t „ = nv u /li is the 
minimum excitation frequency for mode v. 

Note that if we consider the open boundary conditions at the sample ends, x — 0, and, 
x = L, (instead of periodic boundary conditions that are commonly used) then operators 
7] s in Eq.(l) are the same for electrons going right and left. In this case, the electron field 
operator, ip s (x) = ip s +(x) + ip s -{x), vanishes at impurity positions, x = Xi, and expressions 
for the phase fields between the impurity sites turn out to be similar to those found for 
bounded ID Luttinger liquid in Refs.[22, 23, 24], the main difference being the summation 
over j < i that insure proper commutation relations between the electron operators related to 
different segments between impurities. Thus the system breaks up into a set of independent 
segments described as bounded Luttinger liquid with discrete spectrum. 

Now we consider the role of interchain hopping adding to (2) the hopping Hamiltonian 



where indices n and m denoting the chain numbers related to the nearest neighbors are 
added. 

Arguments by Schulz [8] on instability of the Luttinger liquid in the presence of the inter- 
chain hopping were based on calculations of temperature dependence of the thermodynamic 
potential at low temperatures. So we calculate contribution of the interchain hopping to the 
thermodynamic potential per unit volume given by the standard expression [25] 



where V is the volume, T T stands for imaginary time ordering, and (• • •) means thermody- 
namic averaging over the unperturbed state. 

At temperatures T ^> u>i t „, the discreteness of the excitation spectrum can be neglected, 
hence, according to Refs. [6, 7, 8, 9, 10], interchain hopping is expected to give significant 




A-.n) + sin(A, 



-r,m 



A^ n + 2irk F x)\ , (5) 
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contributions destroying the Luttinger liquid. We examine the opposite limit, T <C u>i,u, 
which does not exist in pure infinite Luttinger liquid. 

Consider first the second order correction in t±. The leading contribution to (S) in Eq.(6) 
is given by items in which the term related to a given chain contains contributions from the 
electrons with the same chirality, r, only, 

J2 ^ - J did2(T r e i[Ar ' m(1)_A -> (1)] e~ i[Ar ^^^ (7) 

m,r,r' 

where 1 = {x 1 ,Ti} and 2 = {x 2 ,r 2 }. Other items in which the terms related to the same 
chain contain contributions from electrons moving both left and right give small contribution, 
and we do not discuss them in details. 

Then we use Eq.(4) in (1) and calculate average in (7) using the relation 

Neglecting small corrections oc exp [— due to Planck's distribution functions, we find for 
the average in the exponent 



(T T A r (l)A r (2)) = \\n 



(cosh z — cosy + y Kv Kv ) 



+ % tan 1 



sin y. 



(9) 



(cosh^ — cosy-Y Kv+Kv ^ \_e z — cosy_ 

where y± = z = ^+^-^\) ( chain indiceg are dropped for brevity here). 

In the integrations over 1 and 2, the leading contributions comes from region 1 « 2, i. 
e., |y_| <C 1, z <S 1 where expression (7) reduces to 

t\mL r cos[(r — r')kp{xJ)} dx_dr_ . . 

where e v = 5 V = \{K V + ^-/K v — 2), t_ = \t\ — r 2 |, rr_ = X\ — x 2 , and m is the number 
of the nearest-neighbor chains. Contribution to expression (7) from integration over region 
^ 1, z ^ 1, is small, ~ (a//j) 2<5 , 5 = 5 P + 5 a , because 8 is not too small in the assumed 
case of the not too small interaction (c/. Eq. (3)). 

Additional items in (S) in Eq.(6) in which the terms related to the same chain contain 
contributions from electrons moving both left and right is smaller than that given by Eq. 
(10) by factor ~ (a/l i ) K > ,+K ' T . For reasonable values of K u , this contribution is small and 
can be neglected. 

Similarly, the leading contribution to AQ from higher-order terms in series expansion 
of the exponential in Eq.(6) was found to come from even powers 2n in t± that can be 
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represented as a sum of (2n — 1)!! items like (7) with almost coinciding times and coordinates. 
Therefore, summing up the leading contributions and inserting them into Eq.(6) we can 
calculate the variation of the thermodynamic potential per single chain and per unit length 

t]m r (1 + cos2k F ax) dxdr . . 

= ~ a ^v ' a = J rw((i + rjK v f + x *y+*» ■ 

For moderate repulsion, 8 ~ 1, a ~ 1. In the limit of strong repulsion, K p <C K a ~ 1, a is 
small, a ~ K 2 . 

Thus Af2 is much smaller than the thermodynamic potential of purely ID Luttinger 
liquid, VLq ~ + ^) 

Afi/fio~ (^) 2 , 

and temperature-dependent corrections to Eq.(ll) are determined by small thermally acti- 
vated contributions oc exp [— ^r-]. 

Now we calculate modification of the one-particle Green's function due to the interchain 
hopping. 

G(1,1') = -(^(1M1')5)/(S). (12) 

Again, we consider the low-temperature limit, T <C non-existing in pure infinite system. 
Consider first the second order correction in t± to the Green's function of pure ID system, 

Go(l,l')- 

G 2 (l, 1') = -(T^(l)^(l')^) + (T T i/;(l)j(l'))(S 2 ). (13) 

Calculation is similar to that considered above, (cf. (7-9)). However, in contrast to the 
case of the thermodynamic potential where the leading contribution was given by regions of 
almost coinciding values of times and coordinates, such contributions from two terms in (9) 
cancel each other. So the second-order correction is estimated as 
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G2(1 ' 1 '»~(^) (if^ 1 ' 1 ')' 



Estimation of the fourth order correction in t± gives G 4 ~ (t±l /v F ) 2 (a/ 1) 25 G 2 - Therefore, we 
conclude that at T uji >v the interchain hopping gives small corrections to the one-particle 
Green's function, provided 



'tj\ 2 fa\ 25 ft ± \ 2 /a^ 25 - 2 



where we estimated the cut-off parameter as a ~ 1/kp- 



« 1, (14) 



So far we considered the TL model in which interaction is described by coupling constants 
related to forward- and backscattering. In order to make comparison with experimental 
data we must consider a more realistic Coulomb potential. It is reasonable to assume that 
the long-range part of the interaction is dominated by the Coulomb potential, while the 
backscattering is described by relatively small coupling constant g±. This enables us to 
concentrate on the spin isotropic case and to ignore the possibility of the spin gap. The 
problem of the long-range Coulomb potential on an array of chains was solved in Refs. 
[11, 12, 13]. It was found that interaction of electrons on a given chain is screened by 
the electrons on other chains, and the Coulomb interaction can be described by the TL 
Hamiltonian with coupling constants dependent on wave vector, 

Aire 2 

92 = 94 = s^+e ±q iy (15) 

where s is the lattice period in the direction perpendicular to the chains, and is a back- 
ground dielectric constant for the transverse direction. Coupling constants in spin channel 
remain unaffected. In principle, the coupling constants must be determined by matrix 
elements of Coulomb potential that depend on details of wavefunctions and on chain ar- 
rangement, and must contain an infinite sum over transversal reciprocal lattice vectors. So 
expression (15) is not universal and depends on material. 
Eq. (15) leads to q-dependent velocities 



vf 



K 8e 2 
^ n = —q\\, — = , 1+ 2 , 2 2V K = ^ = Z — ' ( 16 ) 

where k is the inverse Thomas-Fermi screening length. 

We do not perform explicit calculations restricting ourselves to estimations. For the case 
of q-dependent coupling expressions for the thermodynamic potential and Green's functions 
contain various integrals of correlation functions over q^. One can show that the results 
obtained above can be generalized qualitatively to the case of long-range Coulomb interaction 
if we substitute qj_ in Eq. (16) for its characteristic value, q± ~ tt/s. For example, integrals 
for corrections to the thermodynamic potential are dominated by close values of coordinates 
and times, similar to Eq. (10), and coupling parameters should be substituted by their 
averages over q±, 
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5 = -(K p + 1/Kp + K a + 1/K a - 4) ~ - 



(2tt)' \ Kp+ K P/ 



8 



Note that K = 8cx(c/vf), where a is the fine structure constant. Since vf is much smaller 
than the velocity of light, the factor K is large. For vf ~ 2 x 10 7 cm/s, which is typical 
value for transition metal trihalcogenides, K ~ 90. This corresponds to the case of strong 



interaction and leads to quite large values of coupling parameters, 1/K p ~ y ^/e_i_ ~ 3 4 8, 
5 ~ 1/2 4-2. 



Now we discuss conditions for observation of the Luttinger liquid in quasi-lD conductors 
stabilized by impurities. First we discuss the condition for the temperature limiting from 
above the region where the Luttinger liquid can exist. This condition reads, T <C uo^ u . The 
minimal excitation energies, u>i !U , can be estimated as 



where q stands for dimensionless impurity concentration corresponding to number of impu- 
rities per one electron. As Fermi energy in NbSe 3 and TaS 3 is about 1 eV, we obtain that 
ui jP is about 100 K for impurity concentration q ~ 10~ 2 4 10~ 3 . 

Another condition to be fulfilled is smallness of corrections to the Green's function due 
to interchain hopping. According to Eq.(14) the corrections are small provided 



If the interaction is strong enough, 5^1, this condition is not more strict than the condition 
for the limiting temperature discussed above. For lower strength of interaction, 5 < 1, this 
condition reads 



Estimating t± as being of the order of the Peierls transition temperature, T P ~ 100 4 200 
K ~ 0.01e_p, we find that this condition can be fulfilled easily even at 5 = 1/2 for relatively 
small impurity concentration, q 10~ 2 . 

Thus we find that Luttinger liquid can be stabilized by impurities in relatively pure linear- 
chain compounds at rather high temperatures corresponding to experimental observation 
Refs.[15, 16, 17] of the transition from metallic to non-metallic conduction characterized by 
power law dependencies of conductivity and by conductivity resembling the variable-range 
hopping. However, in order to make detailed comparison with the experimental data, calcu- 
lation of the conductivity in a random network made of weakly coupled bounded Luttinger 
liquids is needed. 
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